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ABSTRACT 

Analysis of a 30,000 s X-ray observation of the Abell 3266 galaxy cluster with the ACIS on 
board the Chandra Observatory has produced several new insights into the cluster merger. The 
intracluster medium has a non-monotonically decreasing radial abundance profile. We argue that 
the most plausible origin for the abundance enhancement is unmixed, high abundance subcluster 
gas from the merger. The enrichment consists of two stages: off-center deposition of a higher 
abundance material during a subcluster merger followed by a strong, localized intracluster wind 
that acts to drive out the light elements, producing the observed abundance enhancement. The 
wind is needed to account for both an increase in the heavy element abundance and the lack of 
an enhancement in the gas density. Dynamical evidence for the wind includes: (1) a large scale, 
low surface brightness feature perpendicular to the merger axis that appears to be an asymmetric 
pattern of gas flow to the northwest, away from the center of the main cluster, (2) compressed 
gas in the opposite direction (toward the cluster center), and (3), the hottest regions visible in 
the temperature map coincide with the proposed merger geometry and the resultant gas flow. 
The Chandra data for the central region of the main cluster shows a slightly cooler, filamentary 
region that is centered on the central cD galaxy and is aligned with the merger axis directly 
linking the dynamical state of the cD to the merger. Overall, the high spectral/spatial resolution 
Chandra observations support our earlier hypothesis (Henriksen, Donnelly, & Davis 1999) that 
we are viewing a minor merger in the plane of the sky. 



Subject headings: galaxies: clusters: individual (Abell 3266) 
X-rays 



galaxies: intergalactic medium - galaxies: 



1. Introduction 

Galaxy clusters show a wide range of complex 
phenomenon in the X-ray regime including ra- 
dial abundance gradients, radial temperature gra- 
dients, cooling flows, two-phase gas in the central 
region (Makishima et al. 2001), non-thermal emis- 
sion, and two dimension features such as shock 
fronts. 

Observations with ROSAT, ASCA, and Chan- 
dra have shown that cluster evolution proceeds 
through mergers and that complex physical pro- 
cesses account for the presence or absence of many 
of these observed phenomenon. With the data 
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sets now becoming available, we are in a posi- 
tion to address the physics of cluster evolution. 
ASCA observations of cluster mergers, such as 
Abell 754 (Henriksen & Markevitch 1996) and 
Coma (Watanabe et al. 1999), offer a coarse de- 
scription of the pattern of shocked gas that gener- 
ally matches the merger scenarios played out in N- 
body hydrodynamical simulations. However, with 
increased resolution, Newton-XMM and Chandra 
have begun to provide results that go beyond the 
current simulations of cluster evolution and may 
perhaps require additional physical processes. For 
example, the distribution of metals in the intra- 
cluster medium (ICM) is an important diagnos- 
tic for galaxy evolution within clusters and for 
the cluster gas dynamics. Two-dimensional abun- 
dance maps are within the capabilities of both 
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Newton-XMM and Chandra, yet the distribution 
of metals in the intracluster medium has not been 
simulated. 

Abell 3266 has become a well studied cluster 
merger first shown to have evidence for a merger 
(in addition to optical substructure) based on de- 
tection of shocked gas through detailed model- 
ing of ASCA observations (Henrikscn, Donnelly, 
& Davis 1999). Mergers arc not always apparent 
in the optical data and clusters may even appear 
to have a "relaxed" morphology, either spatially 
or dynamically. For example, though A3266 has 
spatial substructure, it looks dynamically "quiet" 
in the optical in that it has a velocity distribu- 
tion that is statistically consistent with a single 
Gaussian. On the other hand, the ASCA temper- 
ature map and the Chandra temperature map un- 
ambiguously show high temperature regions that 
may be associated with shocked gas. With a de- 
tailed statistical analysis of the galaxy positions 
and redshifts obtained by Quintana et al. (1996), 
in an earlier paper we found that the cluster is 
composed of a secondary component comprised 
of approximately 30 galaxies (~10% of the total 
number of galaxies with measured redshifts) . The 
velocity dispersion is significantly higher, ^1300 
km s -1 in the center compared to the global value 
of 1000 km s -1 . This was interpreted as evidence 
of a merger in progress. Further evidence of a 
merger comes from the co-alignment of the pri- 
mary and secondary galaxy condensations with 
the elongation in the X-ray contours first noted 
by Mohr, Fabricant, & Geller, (1993). The elon- 
gation extends from the arc minute scale in the 
ROSAT PSPC image to the location of the sec- 
ondary galaxy condensation, 15 ' to the NE, and 
is visible in the ASCA GIS image. This, together 
with the Gaussian velocity distribution, suggests 
a merger in the plane of the sky along this axis 
of elongation. The temperature discontinuity in 
the ASCA temperature map is perpendicular to 
the axis connecting the subclusters suggesting that 
the subcluster to the NE has already traversed 
the core and partially shockheated the ICM. Ap- 
plication of the pre- and post-shock gas temper- 
atures derived from the ASCA temperature map 
to the shock equations give a post-shock gas ve- 
locity of 1400 km sec -1 . Simulations of off center 
cluster mergers occurring into the plane of the sky 
(Takizawa 2000), show a similar temperature mor- 



phology to the ASCA map of A3266 but require 
a larger secondary component relative to the pri- 
mary (1:4) than results from our galaxy compo- 
nent separation (^1:10). Thus, the determination 
of merger geometry is somewhat ambiguous using 
only the 2-dimensional temperature map. How- 
ever, a secondary component this large falling into 
the plane of the sky should appear as a spike in 
the velocity distribution and significantly distort 
the Gaussianity of the velocity dispersion. This is 
not seen in our statistical analysis that shows the 
inner region, outer region, and entire cluster are 
all consistent with a single Gaussian profile with a 
similar heliocentric velocity peak (Henrikscn, Don- 
nelly, & Davis 1999). Additionally, there is a nar- 
row angle tail (NAT) and a wide angle tail (WAT) 
radio source with velocities well within the cluster 
velocity dispersion and located, in 2-dimensions, 
within the shocked gas. They are oriented with 
their tail aligned with the direction of the gas flow. 
Therefore, a major merger into the plane of the sky 
seems less likely than a minor merger in the plane 
of the sky. This point is discussed in more detail 
in light of the new results from Chandra in Section 
4.2 

BeppoSax observations showed that the cluster, 
overall, exhibits a radially decreasing temperature 
profile (De Grandi & Molendi, 1999). Chandra 
provides a more detailed description of the central 
region and gives a temperature map that is not 
subject to the PSF deconvolution techniques used 
with ASCA data (Section 3.1). The BeppoSAX 
data did not show any departures from constant 
metallicity out to 10'. Using Chandra, we are able 
to resolve the abundance within the central resolu- 
tion element of the BeppoSax observations, ~2 ar- 
cmin (Section 4.1), and uncover an enhancement. 
These results support the merger hypothesis in 
Abell 3266 and give a more detailed picture of the 
dynamics of the ICM during a merger. Finally, 
the cD galaxy in Abell 3266 has multiple nuclei. 
Multiple nuclei are found in 28% of the first rank 
galaxies in rich clusters (Hoessel 1980). The Chan- 
dra observations provide unique evidence that the 
formation of the multiple nucleus cD in Abell 3266 
may be linked to the merger (Section 4.3). All pa- 
rameters are quoted with 90% confidence errors. 
A Hubble constant of 60 km s™ 1 Mpc -1 is used 
throughout. 
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2. Observations and Analysis 

Abcll 3266 (also known as Sersic 40-6), located 
at (J2000) 04h30m22.0s, -61d32m03s with rcdshift 
of 0.0594, was observed on July 26, 2000 for 30,000 
seconds. The cluster is centered in chip 1 of the 
ACIS-I detector and the emission spans all four 
chips. The total count rate is 10.56 ±0.02 cts s _1 , 
well below the telemetry saturation limit. Figure 
1 is a color coded intensity map that shows the 
full 22' square image in the 0.3 - 8 keV band. The 
image was obtained using processing routines out- 
lined in the CIAO 2.1 Science Threads. The events 
data were first cleaned for flares. The intensity 
map was then background subtracted and divided 
by the exposure map. The background map was 
created following the CIAO 2.1 Science Thread, 
"ACIS Blank-Sky Files" and uses as its basis a 
blank-sky dataset compiled by Maxim Markevitch. 
The background map and intensity map were nor- 
malized by the exposure times of the data sets. 
The normalization was verified by comparing the 
high-energy tail of the background events with the 
tail in the distribution for the study data, which 
should have little emission at the highest energies. 

This intensity map is quite complex and shows 
several interesting features. Overall, the cluster is 
similarly elongated as in the ROSAT and ASCA 
GIS images in the direction toward the secondary 
optical subcondensation to the NE. To the NW, 
there is a steep drop in intensity over 2-3'. The 
intensity peak is visibly displaced to the NW. The 
ROSAT PSPC emission peak was within 30 arc- 
sec of the cD galaxy but the R OS A T pointing er- 
ror did not allow a clear identification. Figure 2 
shows the X-ray contour map from the ACIS ob- 
servation overlaid on the DSS image. The con- 
tours are 1, 1.5, 2, 2.5, 3, 3.5, and 4xl0~ 4 cts 
cm~ 2 s _1 arcmin~ 2 . It is apparent from this fig- 
ure that the X-ray peak is co-spatial with the cD 
galaxy, (J2000) 04h31ml2.16s, -61d27ml5.3s. 

Spectra were extracted from the events file in 
regions of interest. The spectra are grouped so 
that there are no fewer than 64 events per grouped 
bin. A weighted spectral redistribution matrix 
(RMF) and effective function (ARF) for each 
region was found using the calcrmf and calcarf 
tools, version 1.04. The mean RMF (and ARF) is 
computed from the weighted average of the RMFs 
(ARFs) at the different positions in the region. 



The weights are derived from a coarse-grained im- 
age of the data in the region of interest. 

3. Results 

A thermal model was fit to the spectrum ex- 
tracted from a circular region of 2.5' centered 
at the cluster center, (J2000) 04h31m21.04s, - 
61d27m00.16s, for comparison to previous obser- 
vations. The results are given in Table 1. The 
temperature range, 9.17 - 9.60 keV, is consistent 
with the central region as observed by BeppoSax 
(De Grandi & Molendi 1999). The Chandra cen- 
tral abundance, 0.19 - 0.33 (90% confidence), is 
consistent with the BeppoSax value, 0.22 - 0.34 
(68% confidence). The emission weighted tem- 
perature obtained from joint fitting the ROSAT 
PSPC and ASCA GIS observations taken from a 
circular region of radius ^17 ' is 7.75 - 8.46 keV 
while the abundance is 0.15 - 0.25 Solar (Hcnrik- 
sen, Donnelly & Davis 1999). The combination 
of PSPC and GIS gives a similar energy band to 
the ACIS, similar sensitivity above 2.5 keV, and 
significantly higher sensitivity below 2.5 keV. The 
lower temperature of the large region compared 
to the small central region is indicative of non- 
isothermallity; in particular, it requires cooler gas 
outside of the central region. Radial temperature 
(see Figure 3) and abundance (sec Figure 4) pro- 
files were made centered on the ACIS X-ray peak. 
The abundance and temperature measurements 
corresponding to the annular regions are given in 
Table 2. Cooler gas at large radii is detected as 
expected from the earlier comparisons with the 
PSPC/GIS spectra. The temperature decreases 
steeply at >7' as opposed to smoothly decreasing 
as would be expected from a gas in hydrostatic 
equilibrium. The steep temperature ridge may co- 
incide with the edge of the shock front. The Chan- 
dra radial temperature profile in the central region 
is in general agreement with the BeppoSax obser- 
vation, however, the ACIS analysis allows smaller 
radial bins and shows a shallow increase in temper- 
ature within the central 3'. The central 1' radius 
circle is centered on the brightness peak and likely 
has a contribution from the central cD. Cooler gas, 
^1-2 keV, typical of a cD galaxy and the hot in- 
tracluster gas likely account for the lower observed 
temperature. 

The radial abundance map shown in Figure 4 
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spans a radius of ~10'. The best fit Sulphur abun- 
dance (1.0 +/- 0.7 Solar) is consistently ~3 times 
higher than the average abundance of the ICM 
(0.26 +/-0.07 Solar) though they are consistent 
within the large error of the Sulphur abundance. 
The Sulphur-to-Iron ratio is predicted to be quite 
high for SNe II enrichment compared to that of 
type SNe la (along with Si, Ne, O). However, it 
appears to occur at energies just above the Iridium 
edge of the Chandra mirror around 2.2 keV where 
the effective area shows a sharp discontinuity in 
the ACIS response matrix that may depress the 
continuum near the redshifted sulphur line rais- 
ing the inferred Sulphur abundance. Therefore, 
the metallicity was measured tying all of the ele- 
ments except Sulphur which is left as parameter. 
The radial abundance measurements are inconsis- 
tent with constant abundance. There is an overall 
constant abundance with the second annulus from 
the center enhanced at double the cluster aver- 
age. The possibility that the abundance in the 
inner region is low due to multi-phase gas that is 
poorly modeled with a single component as was 
found for M87 ( Molcndi & Gastaldello 2001) was 
investigated. Refitting the spectrum with a sec- 
ond temperature component does not improve the 
fit. The spectrum and residuals from fitting this 
region with a single temperature component are 
shown in Figure 5. If a significant cool temper- 
ature component were present and not properly 
modeled, there should be significant positive resid- 
uals around FeL from the cool component. This is 
not seen in the spectrum. The emission was also 
modeled in 4 quadrants to look for asymmetries 
in the abundance, temperature, and emission inte- 
gral (Figure 6). These results are given in Table 3. 
The temperature does not exhibit azimuthal vari- 
ations. The abundance has large error bars, 30 - 
50%; variations of >50% are not seen. Only the 
emission integral shows variations of 50% due to 
the visibly asymmetric gas distribution. The col- 
umn density is left as a free parameter and is sig- 
nificantly higher than the measured galactic col- 
umn density, 1.6xl0 20 (Dickey & Lockman 1990). 
We refit the spectrum excluding data below 1.5 
keV to minimize the effect of galactic absorption 
on the spectrum. The best fit temperature and 
abundance are consistent with those measured us- 
ing the data above 0.5 keV. The high measured 
galactic column densities do not seem to bias the 



values obtained for the temperature and the abun- 
dance. We also found, by fitting regions spanning 
the entire ACIS-I array, that the high fit n# is not 
local to the central region of the cluster but spans 
the entire cluster region or CCD array. Thus, the 
high uh values may be due to a poor calibration 
below 1.5 keV for this observation. 

3.1. 2-Dimensional Temperature Map 

The ACIS-I observation spans a large part of 
the central and western side of the cluster as shown 
in the temperature map with X-ray isointensity 
contours overlaid (see Figure 7). The regions are 
numbered with the corresponding temperatures 
given in Table 4. The temperature map is color 
coded to facilitate recognition of broad trends. 

To determine the regions to be analyzed, the 
events data for the chips 0-3 were processed 
through a simple algorithm which generates a set 
of regions with a quasi-uniform number of events 
per region. For a given number of events, N, the 
algorithm tests if more than AN events occur in a 
region. If there are more than this number, then 
the region is split into 4 equal-area regions. The 
regions in Figure 7 were created using N — 4096. 
Only four of the 28 regions created had less than 
4096 events (12, 11, 15, and 7) with region 12 
having the least with 2632 events. Nine regions 
have 4000 to 6000 events, and the remaining 15 
regions are evenly spread among having 6000 to 
15000 events. 

A qualitative comparison with the ASCA tem- 
perature map in Henriksen, Donnelly, & Davis 
(1999) can be done since regions 5,6, and 8 on 
the ASCA map are resolved into smaller regions 
in the A CIS map. The hottest region in the ASCA 
map, 11.0 - 16.7 keV, is region 8 which is a small 
section to the southwest of the X-ray peak. The 
A CIS map shows that the primary contributions 
to ASCA region 8 are two regions, with tempera- 
tures at 11 and 12 keV, color coded pale yellow and 
white, respectively. Region 6 in the ASCA map is 
9 - 10.5 keV and is located to the southeast of the 
X-ray peak. While the ASCA map extends several 
arc minutes below the overlapping region on the 
ACIS, the ACIS shows 4 contributing regions with 
temperatures of 7, 8, 10 and 11, keV. The average 
of these 4 regions, 9 keV, is consistent the temper- 
ature of the corresponding ASCA region. Region 
5 on the ASCA map is slightly cooler, 8.6 - 9.6 and 
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extends to the northwest of the X-ray peak. While 
the ACIS map shows a couple of small hot regions 
with temperatures of 11 and 12 keV, however, the 
dominant contribution is from a large region with 
an average temperature of 8 keV. Interestingly, re- 
gion 3, cooler regions in the ASCA temperature 
map, 7.4 - 8.5 keV, is slightly east of the center 
and extends to the northwest. The ACIS map 
shows that this region is dominated by two re- 
gions at 7.3 and 8.5 keV. Thus, it is apparent that 
the temperature map derived from the ACIS, with 
its better spatial resolution, is consistent with the 
temperature map derived from ASCA. The geom- 
etry of the shock front is similar in both the ACIS 
and ASCA maps and generally runs from north 
to south and is centered to the west of the X-ray 
peak. 

4. Discussion 

4.1. Abundance Gradients and Gas Mix- 
ing 

Cluster abundance gradients, including Perseus 
(Ezawa et al., 2001), Scrsic 159-03 (Kaastra et 
al. 2001), A496 (Dupke & White, 2000), AWM 
7 (Ezawa et al. 1997), and the mean metallic- 
ity profile for cooling flow clusters (De Grandi & 
Molendi, 2001), are usually described as monoton- 
ically decreasing radial profiles. Abundance gradi- 
ents may preferentially occur in clusters with cool- 
ing flows (Allen & Fabian 1998) and the BcppoSax 
survey of De Grandi & Molendi tends to support 
this idea. Cooling flows are generally thought to 
be disrupted by mergers. Gas mixing due to a 
merger would erase the gradient as it disrupts the 
cooling flow so that abundance gradients and cool- 
ing flows should anti-correlate with a merger. 

However, there are several clusters with a non- 
constant abundance (reported in the BeppoSax 
survey) without cooling flows: Abell 119, Abell 
3627, and the Coma cluster. Abell 119 has a very 
similar radial abundance profile to Abell 3266 (as 
reported here) while Abell 3627 and Coma appear 
to have an abundance deficit in the center of an 
otherwise flat profile. The abundance change seen 
in Figure 4 shows that there is a significant (>90% 
confidence) difference between the second annu- 
lus and the third, fourth, and fifth annuli. This 
appears to be an enhancement to an otherwise 
flat abundance profile. In the case of Abell 3266, 



the subcluster appears to have deposited sufficient 
metals to raise the observed abundance to 0.4 So- 
lar in the second annulus. Dumping these metals 
would also increase the gas density in this annulus. 
We used the archival ROSAT PSPC observation of 
13,288 seconds to obtain a radial profile and per- 
formed a beta-model fit to it to look for the density 
enhancement. The beta- model is given by S(r) = 
S c (l + (r/a) 2 )<- 3/3+1 / 2 ). The surface brightness is 
described by S with the central value given by S c . 
The core radius (a) and (3 are the other free param- 
eters in the fit. Our temperature map and radial 
profile indicate that the cluster gas is dominated 
by gas at >6 keV. The PSPC surface brightness 
should be insensitive to temperature changes in 
the gas and accurately map the density changes. 
The image was flatficldcd using the exposure map 
and background was subtracted before extracting 
the radial profile. Figure 8 shows a beta model fit 
out to a radius of 10' (~ 1 Mpc) from the clus- 
ter center. The beta model provides a statistically 
acceptable fit, x 2 ~1) to the data within this re- 
gion, however, fitting data at larger radius gives a 
very poor fit due to the pronounced asymmetry in 
the surface brightness that increases with radius. 
Confidence contours for beta and the core radius 
are shown in Figure 9. The range of beta and 
core radius are typical of clusters observed with 
ROSAT (Vikhlinin et al. 1999). In the annular 
region where the abundance enhancement is ob- 
served with the ACIS, r = 0.9 - 1.8 arc min, no 
significant surface brightness enhancement is seen 
above the beta-model. This indicates that there is 
no significant increase in gas density relative to the 
surrounding annular regions. An increase in abun- 
dance without a change in gas density requires 
that metals are dumped in and retained while the 
bulk of the gas is transported away. 

We suggest that the energy input from shock 
heating would preferentially transport the light 
ions (e.g., H and He) out of the subcluster pene- 
tration region since the heavy ions are much more 
massive and therefore have a much lower veloc- 
ity. This increases the metal abundance since the 
metal number density is measured relative to Hy- 
drogen. Figure 10 shows a drawing of the heating 
that is visible in simulations of the cluster merger 
using the Hydra code (Couchman, Thomas, & 
Pearce 1995) that we found matches our proposed 
merger geometry (in Section 4.2). As the sub- 
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cluster penetrates the ICM of the main cluster 
along the proposed merger path, it drives shocked 
gas toward the cluster center and toward the NW 
away from the center. In the post-shock gas, light 
ions and electrons will diffuse out more quickly 
accounting for the preferential retention of heavy 
ions in the region of subcluster gas deposition. 
In the presence of a weak, chaotic cluster mag- 
netic field, thermal conduction will still be efficient 
(Narayan & Medvedev, 2001) and this formation 
process for the abundance enhancement remains 
plausible. Additional support for this scenario is 
apparent in the temperature map that shows sev- 
eral of the hottest regions located 1-2' from the 
cluster center at the approximate location of the 
shocked gas. Figure 2 also shows the compression 
of gas along the NW face of the central region as 
the shock front moves toward the SE. Figure 1 
which shows "fingers of X-ray emission" reaching 
to the NW, presumably expanding post-shock gas. 

4.2. Merger Geometry and Viewing 

The ACIS temperature map shows much more 
detail than was visible in the ASCA temperature 
map. Figure 7 shows the X-ray surface bright- 
ness contours overlayed on the color-coded tem- 
perature map which provides a spatial compar- 
ison of the relative gas density to gas tempera- 
ture. Five resolved regions with both high tem- 
perature, > 10 keV, and high density generally 
lie to the north, west, and south of the X-ray 
peak. The highest temperature regions, >11.5 
keV, coincide with a steep density gradient that 
traces out the high entropy, post shock gas and 
generally runs north to south. In contrast, the 
cluster peak and the region to the northeast is 
cooler, ~8 keV, and has a lower density. These 
general characteristics are crucial for determining 
the merger geometry. Roetigger and Flores (2000) 
make predictions for the temperature map assum- 
ing a merger 3 Gyr ago, off axis, at an angle of 
45° to the plane of the sky, with a mass ratio 
of 2.5:1. While their predicted temperature map 
(their Figure 3) is qualitatively consistent with the 
ASCA temperature map, a significant difference 
becomes apparent with the better spatial resolu- 
tion of Chandra. Their model predicts that the 
X-ray peak should coincide with the highest tem- 
perature gas. The ACIS map shows this is not 
the case. We propose our Figure 10 as the merger 



geometry that is consistent with the temperature 
map. Another problem with the proposed merger 
model of Roetigger and Flores (in which a subclus- 
ter is falling into the plane of the sky at 45 degrees 
with the mass ratio of 2.5:1) is that there should be 
a significant asymmetry on the high velocity side 
of the velocity dispersion. The proposed line-of- 
sight component to the bulk infall velocity, ^1800 
km s _1 , far exceeds the velocity dispersion, ~1000 
km s _1 , of the entire cluster (Henriksen, Donnelly, 
& Davis 1999). Referring to our previously pub- 
lished velocity histogram, in the region of 19,600 
km s _1 , the heliocentric velocity of the primary 
cluster (17,804 km s" 1 ) plus the predicted infall 
velocity of secondary component (1800 km s" 1 ), 
there are only ^30 galaxies (in velocity space) or 
<10% of the galaxies. This would seem incon- 
sistent with a 2.5:1 mass ratio for the proposed 
merger that would predict closer to 100 galaxies in 
the infalling component. Finally, the ACIS map 
shows cool gas running along the merger axis (see 
Figure 10). The relationship of this gas to the X- 
ray peak and the cD galaxy, discussed in more de- 
tail below, is most naturally explained by a merger 
in the plane of the sky. 

4.3. Formation of Double Nucleus cD 

The X-ray intensity peak is slightly displaced 
to the NW of the visual centroid of the X-ray 
emission. The X-ray peak is co-spatial with the 
cD galaxy, (J2000) 04h31ml2.16s, -61d27ml5.3s. 
The cD has a secondary nucleus with a relative ve- 
locity of ^400 +/- 39 km -1 . This is much larger 
than the stellar velocity dispersion of the cD nu- 
cleus, ^327 +/- 34 km -1 . In addition, there is 
an asymmetric rise in the stellar velocity disper- 
sion of the cD that peaks at 700 km" 1 (Carter 
et al. (1985) suggesting that the cD is tidally 
disturbed by a very massive object, perhaps the 
cluster merger. Heating in the stellar velocity 
dispersion is greater than the relative velocity of 
the nuclei themselves which makes it unlikely that 
the double nucleus interaction is responsible. The 
Chandra data build the case for the scenario that 
the formation of the dumbbell morphology is likely 
due to the cluster merger. Figure 11 is an X-ray 
contour overlay on the color-coded hardness ratio 
of the central 6 x 12.5' region. The coolest region 
runs along the axis of the surface brightness elon- 
gation running SW to NE. This ^2.5' plume has 
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a nonlinear shape with its apex at the cD galaxy 
location. Alignment of the plume with the merger 
axis provides a link between the cD galaxy and the 
merger. Tidal heating in the cD and the relative 
velocity of the nuclei are so large that they can 
only be attributed to the merger. This plume of 
cooler gas could to be stripped material from the 
cD galaxy, or a disrupted cooling flow centered on 
the galaxies. The second nucleus may have been 
captured from the subcluster on its initial pass 
accounting for the pre-equilibrium, high relative 
velocity of the two nuclei. The standard picture 
for the formation of a double nucleus is that it is 
formed by cannibalizing neighboring galaxies. In 
the case of A3266, the preponderance of data (e.g., 
the high relative velocity, extreme stellar heating, 
and most importantly the plume that links the 
cD to the merger) seem to indicate cannibalizing 
galaxies from a merging subcluster. 

5. Conclusions 

The Chandra observations of Abell 3266 pro- 
vide a detailed picture of the cluster merger. Ev- 
idence was given to support the hypothesis that 
the merger is in the plane of the sky and is a rel- 
atively minor one. These observations show that 
the merger is undoubtedly linked to the dynamical 
state of the dumbbell galaxy morphology. An off 
center abundance enhancement was found that we 
suggest was formed by the merger. The formation 
process consists of deposition of higher metallicity 
material during the merger followed by a cluster 
wind that reduces the gas density while prefer- 
entially retaining the metals. Future observation 
with ASTRO-E2 XRS will be especially interest- 
ing for Abell 3266 since they will provide gas veloc- 
ities that will directly test our conclusions about 
the merger geometry and formation of the abun- 
dance profile. 
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Table 1 

Single Region Temperature Fit. 



kT a (kcV) 


Abundance 13 


S Abundance 


nH (10 20 cm 2 ) Normalization 


8.6 ±0.5 


0.26 ±0.07 


1.0 ±0.7 


9.9±0.8 1.22 ± 0.02 



a Emission weighted temperature 
b Abundance (excluding S) relative to solar 

c Sulphur abundance relative to Hydrogen: 1.62 x 10~ 5 , Anders E. & Grevesse 
N. (1989) 

d fa(fep / n e n ff dV in cm" 5 



Table 2 

Parameters for Annular Regions. 



Radius 


Normalization 81 


Area b 


kT c 


Abundance d 


< 0.9' 


0.00231 ±0.00009 


2.74 


8.90 ±0.92 


0.21 ±0.19 


0.9' - 1.8' 


0.00420 ±0.00007 


7.92 


9.45 ±0.62 


0.38 ±0.14 


1.8' - 3.7' 


0.01030 ±0.00012 


31.79 


10.0 ±0.65 


0.18 ±0.09 


3.7' - 7.3' 


0.01162 ±0.00014 


69.64 


9.56 ±0.38 


0.16 ±0.09 


> 7.3' 


0.00765 ±0.00013 


159.67 


7.09 ±0.42 


0.08 ±0.10 



\*«D(ilz)r I n - n HdV in cm 5 
b arcmin 2 

c Emission weighted temperature in keV 
d Abundance (excluding S) relative to solar 
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Table 3 

Parameters for Azimuthal Regions. 



Quadrant 


Direction 


n H a 


kT 


Abundance 


Normalization* 1 


1 


North 


9.7 ±1.8 


8.69 ±0.76 


0.30 ±0.10 


0.0124 ± 0.0004 


2 


South 


10.6 ±2.3 


9.13 ±1.00 


0.21 ±0.13 


0.0074 ± 0.0003 


3 


East 


10.4 ±1.8 


8.31 ±0.76 


0.29 ±0.10 


0.0108 ± 0.0004 


4 


West 


9.9 ±2.4 


9.42 ±0.97 


0.27 ±0.13 


0.0089 ± 0.0004 



a Column density of Hydrogen in 10 20 cm 2 
/ n e n ff dV in cm" 5 
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Table 4 

Temperature and Abundance Map Data 



Cell kT (keV) Abundance 



1 


8.7 


±1.1 


0.00 


±0.50 


2 


15.7 


±9.9 


0.01 


±1.23 


3 


7.2 


±0.7 


0.11 


±0.31 


4 


8.1 


±1.0 


0.24 


±0.18 


5 


8.7 


±2.0 


0.35 


±0.15 


6 


7.3 


±0.9 


0.35 


±0.47 


7 


10.7 


±2.8 


0.11 


±0.50 


8 


8.9 


±1.9 


0.01 


±0.37 


9 


9.3 


±3.0 


0.19 


±0.34 


10 


10.8 


±3.2 


0.37 


±0.16 


11 


15.0 


±9.9 


0.08 


±0.45 


12 


11.2 


±2.4 


0.35 


±0.36 


13 


10.3 


±7.1 


0.67 


±0.56 


14 


9.0 


±1.2 


0.18 


±0.35 


15 


7.5 


±1.4 


0.04 


±0.18 


16 


7.4 


±2.0 


0.00 


±0.50 


17 


8.5 


±1.2 


0.23 


±0.15 


18 


7.2 


±0.8 


0.14 


±0.15 


19 


6.5 


±1.1 


0.00 


±0.50 


20 


7.2 


±1.2 


0.10 


±0.19 


21 


8.5 


±2.1 


0.06 


±0.34 


22 


7.3 


±3.0 


0.01 


±0.27 


23 


11.7 ±10.6 


0.24 


±0.25 


24 


5.9 ±2.6 


0.04 


±0.50 


25 


3.2 ±2.0 


0.01 


±0.52 


26 


3.8 


±1.1 


0.02 


±0.20 


27 


4.7 ±2.0 


0.00 


±0.65 


28 


5.6 ±1.2 


0.01 


±0.31 
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Fig. 1. — Color coded intensity map of back- 
ground subtracted image of the A3266 cluster on 
the ACIS-I array. Intensity peak is off center. 
Elongation is visible toward the NE. Fingers of 
emission point to the NW. 




Fig. 2. — Contour overlay of the central cluster 
emission on chip II on the DSS optical image. 
Brightness peak is centered on the central dumb- 
bell cD galaxy. 
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Fig. 3. — The radial temperature profile shows a 
significant temperature drop at larger radius. The 
central 0-1 arcmin region is likely 2-phase with 
the cooler cD contribution depressing the temper- 
ature. 
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Fig. 4. — Radial Fe abundance map shows a sig- 
nificant reduction from the 1-2 arcmin annulus 
and the outer, 7 - 10.5 arcmin annulus. 




hanne energy (keV 



Fig. 5. — Residuals are shown for a single ther- 
mal component model fit to the spectrum from 
the inner 0.9 arc min. This region has a lower 
abundance than the second region. No residual 
emission from FeL is visible that would indicate a 
cool, unmodeled gas component. 



Fig. 6. — The temperature and abundance for 
these 4 quadrants are given in Table 3. The Fe 
abundance is slightly higher in 3 of the quandrants 
and the temperatue slightly lower in the NE. 
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Fig. 7. — Color coded temperature map with in- 
tensity contours overlayed. The regions identified 
correspond to the temperatures given in Table 4. 
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Fig. 8.— A Beta-model fit to the ROSAT PSPC 
radial profile out to 10 arcmin. No significant de- 
partures from a single beta-model are seen at ra- 
dius less than 10 arcmin. 
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Fig. 11. — Contour overlayed on hardness map of 
the central region shows a long filament of cool gas 
running through the central cD along the merger 
axis and to the SW. 
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Fig. 9. — The beta vs. core radius confidence con- 
tours for Abell 3266 indicate a range of paramters 
that is typical of clusters within ~1 Mpc. 




Expanding 
shock-heated gas 



Fig. 10. — Shockheating of the main cluster ICM 
from the merger accompanies the deposition of 
high abundance subcluster ICM. 
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